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COMBINED QUADRti POLE-DIPOLE MODEL FOR INLET FLOW 
DISTORTION NOISE FROM A SUBSONIC FAN 
by Marvin E. Goldstein, James H. Dittmar, and Thomas F. Geider 

Lewis Research Center 

SUMMARY 

A combined quadrupole- dipole model has been developed for the noise generated by 
inlet flow distortion in a subsonic fan. A formula is derived for the total upstream- 
radiated acoustic power in each tone as a function of the design parameters of the fan 
and the properties of the inlet flow distortion. Numerical results are obtained for val- 
ues of the parameters corresponding to various quiet fans tested at the Lewis Research 
Center. The analysis is compared with noise measurements taken on a 51- centimeter 
(20- in. ) diameter research fan as well as with those taken on a number of full-scale fan 
stages. Fairly good agreement was obtained. It should therefore be possible to use 
this model to study the noise-reduction potential of the various fan design parameters. 


INTRODUCTION 

The noise generated by single-stage subsonic fans consists of a number of discrete 
tones, predominantly at the blade passing frequency (BPF), superimposed on a broad- 
band spectrum (which is distributed over a wide range of frequencies). The broadband 
noise is probably generated by such factors as interaction of the blades with patches of 
inlet turbulence and irregular vortex shedding by the fan blades. The discrete tones are 
believed to be caused by periodic spatial nonuniformities (stationary flow distortions) 
entering the fan and by the passage of its wakes across the stator. However, if the 
blade passing frequency is large compared with a characteristic turbulence frequency, 
the sound generated by the inlet turbulence interaction will be concentrated around the 
harmonics of the blade passing frequency and the spectrum will appear to contain tones 
of finite width (refs. 1 and 2). Thus, the inlet turbulence interaction may also be a 
source of pure tones. 



On the other hand, it is possible to adjust the ratio of the number of rotor blades to 
the number of stator blades so that the fundamental blade passing frequency tone gener- 
ated by the wake interaction mechanism will not propagate (ref. 3). Many single-stage 
fans have been designed to take advantage of this theoretically predicted "cutoff. " How- 
ever, the results have been disappointing since these fans still exhibit significant funda- 
mental blade passing frequency tones. 

A possible explanation for this apparent disagreement with theory is that the tone is 
being generated by periodic spatial nonuniformities entering the fan. Inlet flow distor- 
tion can result from crossflows, wakes from upstream bends or obstructions in or ahead 
of the duct inlet, and streamwise vortices from nearby surfaces which are sucked into 
the fan. 

When a narrow-band filter is centered (ref. 4) around a BPF tone produced by a fan 
with no rotor- stator interaction, there are significant amplitude and phase fluctuations 
in the output. This fluttery signal is quite different from the steady and piercing tone 
produced by a rotor-stator interaction. It has been argued by Mani (ref. 1) that the 
"flutteriness" suggests that the tone is produced by inlet turbulence. However, the 
flutteriness in the tone could just as well be produced by an unsteady inlet flow distortion. 

It is shown in reference 2 that the tones produced by inlet turbulence will be much 
broader than the experimentally observed tones unless the turbulent eddies are extreme- 
ly large. It would be helpful to measure the frequency of the unsteady flow distortion 
and compare it with the observed width of the BPF tones. 

Additional evidence against turbulence being the source of pure tones is suggested 
by experiments in which the casing boundary layer (wherein the inlet turbulence enter- 
ing the fan is most intense) was removed (ref. 5). These experiments indicate a much 
larger effect of the turbulence on the broadband spectrum than on the pure tone. 

The nonuniform flow entering the fan can produce noise through either a quadrupole 
interaction or a dipole interaction. The quadrupole mechanism is the result of the fluc- 
tuating Reynolds stresses in the volume exterior to the blades produced by the interac- 
tion of the nonuniform inflow with the potential- flow field of the rotor. The dipole noise 
is the result of the fluctuating lift on the fan blades resulting from fluctuations in the 
angle of attack caused by the passage of the blades through the nonuniform flow. 

Until fairly recently, all analyses of fan noise considered only the dipole mechan- 
ism - the argument that the quadrupoles are less efficient emitters than the dipoles 
usually being invoked. However, such arguments are valid only when the tip Mach num- 
ber is small, and "how small" cannot be determined until after the analysis has been 
made. 

The analysis of the flow distortion noise generated through the dipole mechanism is 
similar to the analysis of the interaction of a rotor with the wakes of an upstream stator. 
This problem has been studied by a number of investigators (see e. g. , refs. 3, 6, and 
7), The possibility that inlet flow distortion noise could be generated by a quadrupole 
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mechanism was first pointed out by Ffowcs- Williams and Hawkings (ref. 8). Morfey 
(ref. 9) carried out an analysis which indicated that the quadrupole noise could dominate 
over the dipole noise at Mach numbers as low as 1/2. 

In this report, a combined quadrupole -dipole model is developed for predicting inlet 
flow distortion noise in a high- subsonic- speed fan. The model accounts not only for the 
individual sources but also for their mutual interactions (which could actually result in 
cancellation). It is believed that this is the first combined model that has been devel- 
oped which accounts for the interaction between these two mechanisms. The quadrupole 
part of the model is a generalization of the one given by Morfey (ref. 9) to include the 
effects of higher order modes (both radial and circumferential). A much more detailed 
model of the source is also given. Thus, in the present analysis, all components of the 
quadrupole which contribute to the sound field are included. The rotor potential- flow 
field is replaced by a line of discrete vortices. 

The analysis leads to a formula for the power emitted from the front of the fan in 
each tone as a function of the design parameters of the fan and the properties of the inlet 
flow distortion. Numerical results have been calculated for a range of parameters cor- 
responding to the various quiet fans tested at the Lewis Research Center. The results 
are compared with experimentally measured tones both from a model fan rotor tested in 
an indoor compressor aerodynamic performance rig and from full-scale fan stages test- 
ed on an outdoor rig. 


ANALYSIS 


Basic Equations of Sound Emission for Fan Stage 

It is shown in section 4. 3. 1 of reference 10 that the density fluctuation p'(x , t) at 
the point x and the time t due to the sound emission from a fan in a duct containing a 
uniform axial flow (with velocity U) is governed by the equation 


P’(x, t) = 


=±r[ 

c 0 J-T Jv(r) 


9 2 G 
0y. dy. 


Tt. dy dT + 


_L f T f 80 

c o J-t Jem 3y i 


f. dS(y ) dr 


(1) 


where T is some large time interval (which will eventually be put equal to °° ), Cq is 
the speed of sound, r is the time (associated with the source), and fj is the force ex- 
erted on the fluid by the solid boundaries . The function 
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t|j ( y ? t) = pvfvj + 6..(p> - ej*p’) - e.. 


(2) 


is Lighthill's stress tensor based on the relative velocity 


v i 


= v i 


U8 


li 


(3) 


instead of the actual velocity v^ In this expression p ! is the pressure (measured 
above the ambient pressure Pq), and e^ is the viscous stress tensor. In equation (3) 
it is assumed that the mean flow is in the x^-direction. The surface S|(t) over which 
the second integral is carried out is the surface of the fan blades, and the volume v{t) 
over which the first integral is carried out is the entire region of the duct external to 
the blades. 

The uniformly moving medium, outgoing- wave Green's function G is given by^ 


G(y- t|£.«) 

r. 

m ? n 


m,n 



where k = oj/cq, M = U/Cq, 0 = |1 - M^, and denotes a doubly infinite set of 

eigenfunctions with eigenvalues k determined by solving the Helmholtz equation 

111 « 11 



$ n + K* „ = 0 

ixi y n xxi y xi xxx y n 


(5) 


on the cross-sectional area A of the duct subject to appropriate homogeneous boundary 
conditions on the duct surface. Thus, when the surface of the duct is rigid, we must 
require that 


1 


The overbar denotes the complex conjugate. 
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( 6 ) 


— ^ = 0 on surface S 

3n 

where 3/3n denotes the normal derivative to the duct surface S. The remaining quan- 
tities in equation (4) are 


r 


m,n 



2 d 72 % 


(?) 


and 


V.= V k 2 -^m,n 

where, in order to ensure that G represents the outgoing-wave Green's function, we 
must choose the branch cut of the square root to be as shown in figure 1. 

When the observation point is far enough from the source (i.e. , in the acoustic re- 
gime), the pressure and density fluctuations will be small and we can use the linear 
approximation 


P’ = cjp* (8) 

The first term in equation (1) is a volume quadrupole source adjusted to account for re- 
flections at the duct wall. The second term is the usual dipole source for fan noise. 


Rectangular Duct Model 

In order to model a real fan, the most appropriate cross-sectional shape of the duct 
is an annulus, such as that shown In figure 2(a). The eigenfunctions S m>n will then be 
combinations of Bessel functions. The analysis can be considerably simplified, how- 
ever, by assuming that the duct is "unrolled” into the rectangular strip shown in fig- 
ure 2(b). (The larger the number of blades, the more closely the rotor approximates a 
two-dimensional disturbance pattern with subsonic phase speed. ) The width 6 of the 
rectangular duct is equal to 2n times the mean radius R of the annular duct, and its 
height b is equal to the outer radius minus the inner radius of the annular duct. In this 
case, we must require that <3> m and its normal derivative take on the same values on 
the surface S a as they do at the corresponding points on the surface (periodicity 
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where the upper sign in the exponent corresponds to upstream propagation (x^ < y^) and 
the lower sign corresponds to downstream propagation (x^ > y^). 


Transformation to Rotor Coordinates 

In order to make the limits of integration in equation (1) independent of time, we 
introduce the rf coordinate system by 

^ = 7i + 5 i2 U t T (12) 

where U + denotes the velocity of the blade row and is taken as positive in the direction 
l t 

shown in figure 3. But since for each value of (ypj^)? fj and T|j take on the same 

values at yg = -6/2 as they do at yg = +6/2, they can be extended to continuous peri- 
odic functions from -°° to +°° . And since the Green’s function is also periodic, the 
integration with respect to ^ can be translated by an arbitrary amount. Thus, the 
integration in the first integral of equation (1) can be carried out over any region Vq 
(of fixed shape which coincides with u(r) at some definite instant of time, and the second 
integral can be carried out over the surface S® of a fixed set of blades. Equation (1) 
then becomes 


P*(x, t) 


CJA* 


(rf , t) d rf dr + 


J C£“' 


(if, r) dS(ff ) dr (13) 


where the linear approximation (8) has been used. 


Derivation of Spectral Equations 

Rather than deal with the pressure directly, it is convenient to deal with its spec- 
trum P(co) defined by the generalized Fourier transform 



(14) 


It therefore follows upon inserting equation (12) into (11) and using the result in equa- 
tion (13) that 
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PM = PqM + p d («) 


(14a) 


where 


CO oo 

Pq(w) = ^ ^ n^^m, n ^ x * ^ 

m=-°° n=0 


(15a) 


oo oo 

PdM = ^ ^ D m,n^ w ^m,n^ x ’ ^ 

m=-°° n=0 


(15b) 


n M = f — — e 

’ i h ^ 


i( Q! 1^7 1 +Qf 2 T ^2 ) 


COS 


D* (w) = 
m s n v/ 


/ 


•j i(aj, 1+ ajr, 2 ) 


L^i 


where 


cos 


(Ol^g) 


l 

r- 


i(w-2mn - U^/6)T 


t! dr drf 


(16a) 


(co-2mffU^./6)T 


± 1 




a l = 


( Mk±k m,n) 


a| = 2mnr/8 
Q!g = nw/b 


f. dr dS(if) 


(16b) 


(17) 


COSj - 


/njTX, 


C* (x , U)) = -i 

m ’ n 2n-6b (1 + 6 n )k 

v n, O' m,n 


exp<-i 




< Mk±k m,n )x l + 


2mffx f 


( 18 ) 


and the plus sign corresponds to upstream propagation while the minus sign corresponds 
to downstream propagation. We shall first treat the quadrupole contribution Pq to the 
pressure spectrum and then the dipole contribution P^. 
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Quadrupole Sound 


Quadrupoie source model . - In order to use these results to predict the sound emis- 
sion, it is necessary to determine the quadrupole source term tT. . As in any aeroa- 
coustic problem, this is extremely difficult to do exactly since it requires a knowledge 
of the complete flow including the acoustic field. It is, therefore, usual to develop an 
approximate model for the source term. As a first step, it is usual in problems which 
do not involve combustion to replace LighthilTs stress tensor t! by the Reynolds 
stress. Thus, we assume 


T ij =p O v H (19) 

We shall further assume that the relative velocity vj can be decomposed into the 
sum of a solenoidal velocity u. associated with the inlet flow distortion and a velocity 
w. associated with the phase-locked rotor pressure field. Then w. is independent of 
time in the rf- coordinate system, and u. is independent of time in the y - coordinate 
system. Hence, 


v?(rf,t) =u.(y) + WjCrf) 


( 20 ) 


where u. and w. denote velocities relative to a coordinate system moving uniformly 
in the y^- direction even though they are expressed as functions of the y- and rf - 
coordinate systems, respectively. 

Representation of distortion velocity . - Since u must be independent of time and 
since it must certainly be periodic in the yg-direction, it can be represented by the 
Fourier series 


u_ 

U 



v ns ys - 
+ ^ B q, s^ cos k + k ^q, s 



2?7iqy /6 
e i 


( 21 ) 


where A. , B„ _, and CL _ are complex constants and i, ], and £ are unit vectors 
in the y^-, yg-, and y^- directions, respectively. In order that u satisfy the solenoidal 
condition 


V • u = 0 


we must require that the coefficients 


£L _ and C„ _ satisfy the condition 

q 3 b H? » 


(22) 
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( 23 ) 




Notice that the normal component of each term in equation (21) automatically vanishes 
on the solid walls at yg = 0 and yg = b. Hence, no further restrictions need be im- 
posed on the coefficients A„ _ and B„ _. Since the problem is linear, we need only 
calculate the sound generated by a single harmonic component 




B q,s sin 


Vsy 3 Vj 27riqy 2 /6 


and the solution to the complete problem can be obtained by summing the results over 
q and s. 

Rotor velocity field model . - We shall now develop an approximate model for the 
phase-locked rotor velocity field Wj. Most fans involve fairly small camber (especially 
at the tip, which is the critical region for producing sound), and thus linearized theory 
can be used. But since the quadrupole source is of interest at high subsonic Mach num- 
bers, we must take into account the effects of compressibility. It is convenient to in- 
troduce the coordinate system which is alined with the blades as shown in figure 4. It 
is well known that the potential- flow field about such a blade row can be represented by 
a uniform velocity plus a distribution of line vortices along the blade chords. In order 
to simplify the analyses, we replace this vortex distribution by a discrete set of concen- 
trated vortices lying at the midchord points. Let and V 2 denote the velocity com- 
ponents along the Xj- and X 2 - directions, respectively (relative to the blades). Then 
according to linearized theory the velocity induced by a single compressible vortex of 
strength r q at the origin is given by (ref. 11) 


Vl-i 


P r J 2ui Z 


where 


^r= A 1 


Z — X| + i/3 r X 2 
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Hence, summing up this result over all blades shows that the concentrated vorticity ap- 
proximation to the velocity field about the cascade is 



(27) 


where Vq is a uniform (in general) complex velocity which will be chosen to make the 
flow far upstream of the cascade equal to U r and 

A o = A ^r cos y - 1 sin y) ( 28 ) 

where A is the interblade spacing and y is called the stagger angle. But upon using 
the relation 


oo 

- coth Z 

Z - inn- 

n=-°° 


equation (27) becomes 



And since 


we must put 


Hence, 


lim coth — = - 1 
x r -eo a q 


lo 

^r 


r 

2A 
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(29) 




+ coth — 
A, 


It can be seen from figure 4 that the Xp X 2 coordinate system is related to the rotating 
rjpTjg duct-alined coordinates by 

Xf = cos y + r ?2 sin y 

X 2 = sin y + f ?2 cos 

Hence, it follows from equations (26) and (28) that 

z 77 ., (cos y - i 0 sin y) rj 2 

— + i — 

Aq A(/3 r cos y - i sin y) A 

_ + i(l - $1) sin y cos y] r ? 2 

* /n2 2 *2 \ A 

A(p r cos y+ sin yl 

But upon using the velocity triangle in figure 4, this can be written as 



where 


z = + ii 7 2 


D * 


/3 r + iMM t 


/3 


2 


and 


M 


t 


U, 


c 0 


(30) 


(31) 

(32) 


(33) 
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is the Mach number of the blade row. Now the velocity components w x and w 2 are 
in the 77 j- and ^-directions and are measured relative to a reference frame moving with 
a uniform velocity U in the ^-direction, while the velocities V x and V 2 are in the 
X j- and X 2 -directions and are measured in a reference frame moving with the cascade 
velocity U t . Hence, it follows from figure 4 that 


Vj = (wj + U) cos y + (w 2 + U t ) sin y 


and 


V 2 = (w 2 + U t ) cos y - (w^ + U) sin y 


= w 2 cos y - w^ sin y 


Then 


(fi r cos y - i sin y) 


V x - if 




- u_ 


w x /3 2 - iw 2 D 


Hence, multiplying equation (29) through by (fi r cos y - i sin y) and using equations (28) 
and (30) shows that 


fi^(w 1 - iw„D) = — ^ (l + coth — \ 
12 2iA V A/ 


(34) 


But since 


r 

lim j3 2 (w., -iw 0 D) = — — fir 
Vl ~+°o V 1 2 > iA 


equating the imaginary parts of this expression shows that Aw 2 , the change in w 2 
across the blade row, is related to r q by 



■¥ 


where 6 = -Aw 2 /U t is called the work coefficient of the fan. Hence, equation (34) can be 
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written as 


2 _ ie/3 u , v 

P (w^ - iWgD) = — - — + coth-2— J 


i^ r u t e nz / A 

2 sinh 15 


By using the geometric series 


— — = ^ Z p for |Z| < 1 

1-2 p=0 


(with Z — exp (±2»rZ/A)), this becomes 


- 10^9 £ ® 2 ’' PZ/4 « Vi < 0 

P =1 


0 (Wj - iw 2 D) = < 


iu t )8 r e e -27rpz/A if Vi> Q 

p=o 


Equating real and imaginary parts now shows that 


Wj/3 2 - w 2 MM t = t 0/3r ^ ( e 2n-pz/A _ g 2ffpz/A ) for ^ < Q 


P=1 


ils 


w. 


= — ^ (e 2l[>z/a + e 2 ’P z / A ) lor ,j < 0 


P=1 


(35) 


or adding the results 


-iU,0 


W-| = 


L £ (De 2lT P z ^ A - De 2 'P z / a ) for ^ < 0 


P=1 


(38) 
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The noise generated by inlet flow distortion probably 1ms its largest effect on the sound 
field passing through the fan inlet. We shall therefore restrict our attention to the 
upstream -propagating waves. We shall also suppose that only the front half of the cas- 
cade potential- flow field is effective in interacting with the inlet distortion and generat- 
ing sound which propagates upstream. Hence, we set 

W 1 =w 2 - 0 for tj^ > 0 (37) 


Evaluation of quadrupole source term , - The results of the previous two subsections 
can be used to evaluate the source term n in equation (15a). Thus, using equa- 
tion (12) to eliminate y. in equation (24) and inserting the result in equations (19) and 
(20) shows that equation (16a) becomes, upon carrying out the integration with respect 
to r)g (since only the cross terms w.Uj can contribute to the sound field - the other 
terms representing steady flows in a inertial reference frames), 


Q^n = -bU( 1 + a S;0 )5 S;n Po< E m 


n <5 

exp 1 

J 


cfyl + y (q + m)r} 2 l w a drj 2 dr\ 1 


x 



exp<i 


r 

2 7TUj. 



a; (m + q) 

T> 

l 

U 

J 


> dr 


(38) 


where (since Wg = 0) the repeated Greek index a takes on only the values 1 and 2 and 


E m,n«* a i A . 


+ “2 B q,s + 


a 3 2 1 b 
s 6 9> s 


Hence, it follows from equation (17) that when n = s 


E (k)=MHUL2) B (Mk + k m,s } A 
m,s w 6 q,s + 2 q, s 


P 


(39) 


Carrying out the integration with respect to r and taking the limit as T — now 
shows, since lim / e iT ^ dr - 2tt 5(|3), that 

T->oo J-T 
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Q™ n = “ b (l + 5 s, 0 )U6 s, n 2n 5 ( w " — + ^ 1 P(\ a 


0 r6 


expji a^j} 1 + ^ (q + m)r] 2 |w a dn 2 dr j l 


Inserting equations (17), (35), and (36) into this result and carrying out the integrations 
with respect to and 772 show that 


b5p 0 U t e(l + 5 Sj0 )6 s?n U,6 co - — (m + q) E 


..+ 2nm 


icet - — (m + q)D 

1 5 


if “JJU-1,-2,-3, 
B 


Q+ m n = ) 0 if HL±_2 4. ±1, ±2, 
’ B 


-b5p o U t0 (l + 6 fl> 0 )6 Sj n U,5(co - 4 1 (m + *>) E m, s 


iDa!j - 


iat - — (m + q)D 
1 6 


if = 1 , 2 , 3, . 

B 


where the overbar denotes the complex conjugate. Substituting equation (40) into equa- 
tion (15a) and changing the index of summation from m to m + q now shows that 


V“ )=M W e < 1 + 6 s,0>" 


/ 2?ru t \ . 

5 k B - q , s (*) 


p=- °° 

p*0 


x E R /2^M t p\ 

p,q,s pB-q, s l A j 
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p,q>s 


H p,q, S for p > 0 


H p,q,s for p < 0 


p>q,s 


_( P 'l)^' iD(MM ‘ P + K P.q,s> 


and (see eq. (10)) 


q,s 2jr pB-q, s 


’2nM^ 


W- A-af +/M \ 2 

1 \ B/ \2b/ 


The pressure spectrum for an arbitrary flow distortion can now be obtained by summing 
equation (41) over all s and q to obtain 


V W) - '“P 0 U t C9 X (1 + « g> „) Y sL - 2 J^Pjc; B _ q s (x) 

p#0 


x E tj /2^rM£\ 

p,q,s pB-q, s l A j K 


This is the formula for the quadrupole source contribution to the pressure spectrum. 
The principal assumptions on which it is based are discussed in the section Quadrupole 
source model. 


Dipole Sound 


We shall now calculate the contribution of the dipole source to the pressure field. 
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Evaluation of dipole source term for periodic forces . - It is usual practice when 
dealing with dipole sound from fans and propellers to neglect the variation in the retard- 
ed time ajT) j between the surface of the blades and the rotational plane of the cascade 
(which we take to be the 77^ = 0 plane). It is also reasonable to neglect the blade forces 
in the (radial) 7)3- direction. Then it is shown in reference 11 that the integration in 
equation (16b) over the blade surface can be reduced to an integral over the rotational 
plane of the cascade to obtain (after carrying out the differentiation with respect to rjj) 


„(«) = iaj 


m,n 



ia &2 

e 



i(o)-2mfflL./6)T^, 
e f a dr dr? 2 


dr) 3 


(45) 


where, as before, the Greek index a only takes on the values 1 and 2 and ^(77^, 773) 
is the a** 1 component of the net blade force per unit projected area (on the rotational 
plane) acting on the cascade at the point 772,7)3 of the rotational plane. 

Since it is assumed that the inlet flow distortion is steady, the unsteady blade forces 
due to the passage of the blades through this distortion must certainly be periodic, with 
frequency equal to the shaft rotational frequency 

2jrU t 

S2 = 5 (46) 

6 

Hence, we can expand these forces in a Fourier series to obtain 

% - T) r >2^3> e ' 1,nT W 

q=-oo 

where 

r&Tr/Cl) n 

F« = ~ f e iqfiT f dr (48) 

* 2 n J 0 ct 

Upon inserting this into equation (45), carrying out the integration with respect to r, 
and taking the limit as T — °° , we get 


m 


n (co)=ia G 2u 


£ 6 (w - (m 


q=- c 



iQ! 2 7? 2 

e cos (Q!3T? 3 )Fq dr) 2 dr) 3 


(49) 
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Inserting this result into equation (15b) and shifting the index of summation from m to 
p = m + q and using equation (17) show that 


oo oo 


P D M=2ri£ £ 6(a,-pn)C< 0 ) qjn (x) 

n=0 p,q=-°° 


i2ff(p-q)rj 2 /5 



x< 


Ml: q) f 2 + — 

5 « p 2 


2jrMM,. 


P ±k p-q,n(°P) 


F J dn 2 dr? 3 (50) 


where 


c (0) 

P,q,n 


(x) = Cp_ q? n (x, ftp) 


(51) 


Case of identical blades . - When the fan consists of B identical blades, the un- 
steady force acting on each blade must be the same as on any other blade whenever it 

passes through the same position of the inlet disturbance. Hence, the force f can be 

0 u 
expressed in terms of the individual blade force f by 


B 


VX f !^ 3 ’ ^2 + | < s - H - !>) 

S=1 


Inserting this into equation (48) and shifting the variable of integration implies that 

F q = Yj e' 2, ' 1(S ' 1>q/B <4(13’ 42 + 5 - d) (52) 

S=1 


where 



^iqftr f 0 


f a^2 ? ^3’ T ) dT 


(53) 


A 

is simply the q Fourier coefficient of the force f^ acting on an individual blade. 

Substituting equation (52) into equation (50), shifting the variable of integration from 
rj 2 to tj 2 + (6/B) (s - 1), and noting that 
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show that 


B 


-2irip(s-l)/B 


s=l 


fiB for p = 0, ±B, ±2B, . . . 
|o for p + 0, ±B, ±2B, . . . 


P D (w) = 2ffiB ^ 
n=0 


z 


p,q=-« 


6(w - pBO)C 


( 0 ) 

pB, q, n 


(x) 



e i2jr(pB-q)r] 2 /6 


COS 



2ff(pB - q) F 0 
6 2 ’ q 


1 

2n\MMj.p 

±k 

,3 2 

A ] 


pB-q,n 


(fipB) 



where we now need only carry out the integration over the projected area Aq of a single 
blade. 

Dipole source model . - We now encounter the same difficulty as occurred with the 
quadrupole term. Namely, the blade force f^, which must be known in order to calcu- 
late the sound field, is, in general, a complicated function whose exact determination 
requires a knowledge of the complete flow, including the sound field. We shall again 
proceed by developing an approximate model for this term. Although it is now possible 

to obtain a fairly sophisticated model for f®, we shall adopt the relatively simple model 

u 0 

that is conventionally used in the literature. Thus, we suppose that f” is concentrated 
at the centerline of the blade (which we can suppose lies at ^ = 0). We shall also sup- 
pose that the viscous contributions can be neglected. Then f® can be expressed in 

terms of the fluctuating lift force per unit span l acting on the blade and the stagger 

2 

angle y defined in figure 4 by 


^2 = cos y 
f 1 = 6(ij 2 )£ sin y 

Inserting the result into equation (52) now shows that 


2 Re call that 
ed on the blades. 


denotes a force exerted on the fluid while 


denotes a force exert- 


20 



> 


( 55 ) 


F 2,q = ' 6( ’ , 2 )L q COS ’' ! 

tVi 

where denotes the q Fourier coefficient of the lift force per unit span acting on 
an individual blade. 

Up to this point the assumptions have been fairly mild. However, we must now de- 
termine Lq from a knowledge of the coefficients of the flow distortion entering the fan. 
We shall restrict our attention to the case where the flow distortion is only in the axial 
velocity so that B„ = C = 0. Since the quadrupole term dominates the dipole at 

4? ® 4? ® 

sufficiently high Mach numbers, we shall not take into account compressibility in calcu- 
lating L . In addition, suppose (1) that each airfoil in the cascade acts individually and 
is uninfluenced by the remaining blades, (2) that it responds locally in a two-dimensional 
manner so that strip theory can be applied, and (3) that the camber is small enough so 
that linearized thin- airfoil theory applies. 

Then, as is shown in reference 10, can be calculated from the Fourier coeffi- 
cients in equation (21) and the Sears' function S(a ) by 

4 

7TST]q 

A a s cos “ ( 5 6) 

q,s b 

s=0 

where c is the airfoil chord, S(ct^) is the Sears' function, and 

cr q *S0£ (57) 

q 2U r 

is the reduced frequency. The Sears' function can be expressed in terms of Bessel func- 
tions. However, it is well known that it can be approximated quite accurately for almost 
all frequencies by the relation 


= -ircf> 0 U r U Bin rS(a^) 


S <V 


exp 



2(1 + 2 kg! )j 



(58) 
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We shall adopt this approximation herein. Bender (ref. 12) has recently shown that 
in the case of a flat plate the correction for finite length in the y^- direction can be ap- 
proximated fairly accurately by multiplying equation (56) by the ratio of the steady- state 
lift slope dC^/dar (the lift slope for infinite-span thin airfoils) to 2tt. Upon the adopting 
of this correction, equation (56) becomes 


Pq dC T V 

— — UU sin *3(a fl ) ) A „ cos 


V 


2 dot 


/ j 

s=0 


TTSTJq 

b 


(59) 


The lift- slope dC^/da is a function primarily of the aspect ratio AR of the blades 


AR = — 


(60) 


tion 


Finally, it is shown in reference 13 that dCjJda can be approximated by the rela- 


dC L 2n AR 

da AR + 2 

\AR + 2/ 


(61) 


Using the relations (shown in fig. 4) 


sin y = — 
U r 


cos y = — 
U r 


substituting equation (59) into equation (55), using the result in equation (54), and carry- 
ing out the integrations show 
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P D (o,) = ri bBAp 0 U t U i 


1 dC L c C 01 


2 da A U 


< 1 + 6 sV 


2 nil. 


w - 


r/s,q 


J P C pB jq ,n^q) 


P=- c 


2r(pB-q) M l/ 2 *MM t 


6 


p 


p±k pB-q,s H 


q? s 


(62) 


Evaluation of Pressure Field Due to Individual Tones 


At this point, we shall restrict our attention to upstream propagation for the case 
where the flow distortion is purely in the axial direction. Then inserting equations (41a) 
and (62) into equation (14a) and using the result in equation (14) show that the pressure 
field only contains frequencies which are multiples of the blade passing frequency and 

in 

that the p harmonic of the pressure fluctuation p’ is 

r 


2bBpgjr 2 U. U /-2irU t \ _ 

P P ~2 ^ ( a P y 2 ^ ^ + 5s > 0^ C PB, q, s^ x )G p, q, s A q, i 


(63) 


where we have used equations (39) and (44) and put 


- n 

®p, q, s ' + K p, q, s^q, s + ‘* S <V f? - + q> s ) 


(64) 


and 


1 /_ \ -1 

i// =- - jM" 1 

2 da \A/ r 


(65) 


Radiated Acoustic Power 


The quantity which is perhaps of most interest is the total acoustic power radi- 

in P 

ated upstream in the p harmonic. This can be calculated by integrating the axial 

in. 

component I (x ) of the sound intensity in the p harmonic over the cross-sectional 

r 
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area of the duct. Thus, 


rd fb _ 

n = lim 2 If I (x) dx 9 dx„ for p=l,2, 

P x Jo JO P * d 


( 88 ) 


It is shown in reference 11 that 


Ip = (1 + M 2 )PpV p + 


“-|Pp | 2 + ^ 0 M|v p | 2 
P 0 C 0 


(67) 


where v is the p harmonic of the acoustic velocity and is related to p' by 
P P 


-api / 2 m, , \ 


Hence, inserting equations (18) and (63) into this relation (with k = 2;rMj.p/A) shows that 

.2 


2bBU, Uff -2irpU t t/A _ 

v„ = e ^ x s , q (1 + 6 s, 0 )C pB, q, s (x )0 p, q, s A q, i 

s,q 


p 2 
c 0 /3 2 


( 68 ) 


where 



MM t P + K p,q,S 

^ +MK p,q,s 


(69) 


Since equation (18) shows that 



dXg dxg = 


6 ,6 

m,m' n, n* 


* 5 „, 0 > 


substituting equations (63) and (68) into equation (87) and using the result in equation (66) 
shows that 

3 — — 

The factor 2 arises because both I and I „ contribute to the acoustic power in 

th P - P 

the p tone. 
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(MM t ) 2 pM t 
p Q c 2 BbA 4 



For convenient reference we also note that 


H p» q, s 



-iD(MM t p + K pqjS ) 


P/3 r 


" lK p, q, s 



© 


p,q 5 s 


-0Hp > q j s(M M tP + Kp 5 q ?s ) + i^S(a q ) 


M(pB - q) ^2 
B 


M t(MMtP + K p, q 


* =• 


271 AR 


2IVL 


AR + 2 


(AR + 4 
IAR + 2 


exp -lcr 


S <V = 


2(1 -5- 27r|cr | ) 


q 1 u 


V 1 + 27 r l CT q|) 


cr 


q 



( 70 ) ,, 


(43) 

(44) 

(71) 

(72) 

(58) 

(73) 
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Equation (71) was obtained by combining equations (42) and (64), equation (72) was ob- 
tained by combining equations (65) and (61), and equation (73) was obtained by combining 
equations (57) and (46), 


Calculation of Distortion Coefficients 

In order to use equation (70) to predict the radiated sound power, it is necessary to 
determine the Fourier coefficients A_ But it follows from equation (21) that these 

t.? ® 

coefficients can be determined from the axial distortion velocities u^ by the relation 


< l’ s b6(l + 6 


-ft 

’’ 0 Jo Jo 


27riqy 2 /5 /VsyA u 1 (y 2 >y 3 ) 


cos 


u 


dy 3 dy 2 (74) 


It is convenient to suppose that the distortion velocity can be expressed as the product 


u. 


U 


= i f(y 2 ) s(y 3 ) 


(75) 


then 


\ s Ia q d s 


where 


a 


q 



-27riqy 9 /5 

e f(y 2 ) dy 2 


and 


d = f cos g(y 3 ) dy 3 

8 b < 1 + 6 s,o)jf 0 \ b J 

When the flow distortion can be represented by N Gaussian profiles (for N = 1, 2, . . .) 
in the circumferential direction (see fig. 5(a) for the case where N = 2), 


26 



f(y 2 ) = exp 


y 2 1 + 2j 


£ <v 0 


for j = 0, 1, . . . ,N - 1 


we find that ' 


VFNa’e N * e'^’) 2 


0 otherwise 


~ n q ^ Aeerf(—L- + Triqo'') for q=0, ±K, ±2N, . . . 
\2Nor’ j 


When the flow distortion is a square wave in the circumferential direction (fig. 5(b)), 


f(y 2 ) = 


1 for 


1 

6 2 


0 for — - A > a" 
5 2 


and we find that 


a =l_^sin2!rqcr” 
H ?rq 


When the flow distortion is uniform in the radial direction 


we find that 


g(y 3 ) = 1 


d s = 6 s,0 


And when the flow distortion is a square wave concentrated around the tip region 


*If a ! gets large, the Gaussians will merge but, of course, still remain periodic. 
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g(y 3 ) = < 


y 3 

for 1 - a < — si 
b 


for 


0 



< 1 - a 


(82) 


we find that 


d 


s 


(-l) s 2 
is(l + 5 s>0 ) 


rs-p 

sin user 


(83) 


RESULTS AND DISCUSSION 
Parameters in Solution 

Equations (70) to (73) can be used to calculate the total sound power radiated up- 

ii. 

stream in the p harmonic of the BPF tone which results from the interaction of the 
rotors with an axial inlet flow distortion. In order to do this, it is necessary to know 
the axial distortion velocity u^ entering the fan. The geometry of the fan is character- 
ized by the parameters c, b, B, and A. They denote the blade chord, the blade height, 
the number of blades, and the interblade spacing. The remaining parameters which ap- 
pear in equation (70) depend on the operating conditions of the fan. They are the cascade 
velocity (tip speed) U^, the throughflow (axial) velocity U, and the work coefficient 6. 
The work coefficient is a measure of the loading on the fan. 


Dipole and Quadrupole Contributions 

The first term in equation (71) represents the quadrupole contribution to the radiat- 
ed sound power, while the second term represents the dipole contribution. These terms 
will not always be additive but can easily be out of phase and produce a sound field which 
is less than that which would be produced by either one alone. It follows from equa- 
tions (71) and (72) that the ratio of the quadrupole term to the dipole term is proportional 
to the ratio of the work coefficient 9 to the cascade solidity c/A. Thus, the quadru- 
pole contribution will tend to dominate for highly loaded, low- solidity blades. The 
dipole contribution is independent of the work coefficient. 
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Numerical Results 


The dimensionless radiated sound power in the first tone due to a localized flow dis- 
tortion has been calculated for values of the fan design parameters covering a range 
corresponding to the various quiet fans tested at the Lewis Research Center. 

It is assumed that the flow distortion consists of a single pulse with a Gaussian 
shape in the circumferential direction and is a square wave in the radial direction (i.e. , 
a combination of eq. (76) with N = 1 and eq. (82)). It has been found, that for narrow 
distortions, the detailed shape of the pulse has little effect on the radiated sound and only 
its size is important. Equations (70), (74), and (75) show that the radiated sound power 
varies as the square of the normalized distortion velocity. 

The plots show the variation in radiated power in the fundamental tone as the rotor 
circumferential Mach number is varied while the ratio of axial to rotor circumferential 
Mach number is kept equal to 1/2. Most quiet fens tested at Lewis correspond to this 
case when operated along a fixed throttle line which passes through the design point. 

When fans are operated in this manner, the work coefficient 9 tends to remain fairly 
constant. In fact, the work coefficient varies little from fan to fan and is roughly equal 
to 0. 4. 

Figure 6 illustrates the effect of changing the number of blades while holding the 
cross-sectional area of the fan constant. It can be seen from the figure that the effect is 
very small (of the order of about 2 dB) over the range considered. 

Figure 7 shows the effect of changing the circumferential size of the distortion. 

The distortion width is 2a'5. It can be shown from the figure that for very small distor- 
tions the sound power varies roughly as M^ to the sixth power, while (at least for this 
case) when the distortion becomes about equal to the blade spacing, the tone power var- 
ies as the velocity to the eighth power. The fans tested at Lewis exhibit anywhere from 
a velocity-to-the-fifth-power relation to a velocity- to- the -eighth- power relation (see 
fig. 20 of ref. 4). As the size of the distortion increases further, the rate of increase 
of the radiated tones with velocity becomes greater. This is due to the fact that the inlet 
distortion contains predominantly lower order harmonics which generate acoustic modes 
which are cut off at the lower Mach numbers. Notice that at each Mach number there is 
a certain distortion size which maximizes the noise. 


Comparison with Experimental Data 

The predictions of the analytical results have been compared with various fans test- 
ed at Lewis. The first comparison is with a 51-centimeter (20- in. ) diameter model fan 
rotor tested in an indoor compressor aerodynamic performance rig. The results of 
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these tests are reported in reference 14. This fan, designated rotor 1 - model 1, is 
shown in figure 8 . The relevant design parameters and operating conditions are listed 
in table I. The latter correspond to the highest subsonic speeds at which the fans were 
tested. The flow distortion entering the fan has not been measured. However, there 
are four support struts ahead of the fan. Downstream of the struts there is a large con- 
traction of the flow to the 51- centimeter (20-in. ) diameter fan size. We therefore as- 
sume that the flow distortion can be represented by four equally spaced Gaussian pro- 
files in the circumferential direction, where the width is equal to the thickness of the 

struts, and that it is uniform in the radial direction. The amplitude of the distortion is 

1 

unknown, but a maximum value of percent of the free-stream axial velocity is not 
unreasonable. The comparison between the experiment and theory is shown in figure 9. 

Since the amplitude of the flow distortion has not been measured, the absolute level of 

1 

these curves has to be regarded as uncertain. Indeed the assumption of percent 
distortion is equivalent to putting the theoretical curves through a single data point. 
However, it can be seen that the theory certainly predicts the correct trends. The trian- 
gular points show the contribution of the quadrupole source to the sound power. It can be 
seen that at the center of the 80- percent- speed curve the quadrupole is the dominant 
contribution, while for the 80-percent-speed curve the largest contribution is from the 
dipole source . The dropoff at the high- Mach- number end of each of the constant- speed 
curves occurs because of the dropoff in the dominant quadrupole source. The dipole 
term tends to decrease monotonieally with decreasing relative Mach number. 

The analysis was also compared with front- quadrant fundamental tone measure- 
ments taken in several subsonic full-scale fans in the installation shown in figure 10. 

The concrete support structure in this installation is believed to be the cause of an inlet 
flow distortion entering the fan. This distortion was measured for the QF-3 fan in ref- 
erence 7. By using this measured flow distortion data as an input to the analysis, the 
fundamental BPF tone is predicted. This prediction is compared with measurements on 
a number of fans operated subsonically on this rig. The relevant design parameter and 
operating conditions are shown in table II. All the fans were nominally 1. 8 meters 
(6 ft) in diameter. 

In carrying out the comparison, it is assumed that the flow distortion is the same 
for all fans as it was for the QF-3 fan (operated at 90 percent speed) for which the 
measurements were taken. The measurements of reference 8 showed that a small re- 
gion of high velocity existed at the bottom center of the inlet. This distortion is modeled 
analytically by a Gaussian profile in the circumferential (yg) direction and a step profile 
in the radial (y^) direction. The measured distortion was found to cover 10° in the cir- 
cumferential direction. Hence, we take the ratio 2a ? of the width of the Gaussian to the 
circumference to be 10 °/ 360 °, or 0. 139. The radial extension of the distortion was 

1 (v 

about 8.9 centimeters (3^ in, ). Hence, the ratio a of the height of the distortion to 
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the blade height is taken to be 3. 5/19, 35 or 0. 18. By using these data, numerical re- 
sults have been calculated for the fans listed in table II. The results are compared 
with the measured levels of the fundamental tone in table 331 . 

It can be seen that good agreement is obtained for the QF-3 fan, on which the dis- 
tortion measurements were made, and for the nearly similar QF-1 fan. However, the 
agreement with fans A and B is not very good. One possible explanation of the discrep- 
ancy is that the distortion entering fans A and B may not be the same as that measured 
on QF-3. Another possible explanation is that the analysis does not take into account the 
back reaction of the rotor potential- flow field on the inlet flow distortion. This back 
reaction could cause the distortion velocity distribution to be deformed before it inter- 
acts with the rotor. As a result the rotor would sense an inlet flow distortion with a 
higher harmonic content, which would then produce a scattered pressure field with more 
higher order modes. Since these modes will not be cut off as readily when the Mach 
number drops, more sound will be radiated at lower Mach numbers. Because fans A 
and B are run at lower Mach numbers than QF-1 and QF-3, the inclusion of this effect 
might improve the agreement with the data from the former fans. 

A third possible explanation is that the outlet guide vanes in fans A and B, which 
are closer to the rotor than they are in QF-1 and QF-3, are producing some kind of 
rotor-stator interaction noise which is not being cut off. Finally, it should be noted 
that much of the energy in the discrete tones may be the result of unsteady flow distor- 
tions which cannot be detected by a velocity traverse made with a slow-response instru- 
ment such as a pitot tube rake. 


CONCLUDING REMARKS 

A combined dipole -quadrupole model (which includes the cross coupling between the 
two mechanisms) for predicting inlet flow distortion noise has been developed. Numer- 
ical results were obtained, and a comparison with the fundamental blade passing fre- 
quency noise from subsonic fan stages tested at the Lewis Research Center was made. 

Very good agreement was obtained with data from a 51-centimeter (20- in, ) model 
fan (tested indoors) by effectively adjusting the level of the theory so that it goes through 
a single point. This adjustment was necessary since the absolute amplitude of the dis- 
tortion (which enters the equation as a multiplicative factor) was unknown. Good agree- 
ment was also obtained with the level of the fundamental tones measured on full-scale 
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fans. In this case, no adjustment was necessary since the distortion was measured. It 
should therefore be possible to use this model to study the noise reduction potential of 
the various design parameters of the fan. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 7, 1974, 

501-04. 
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APPENDIX - SYMBOLS 


AK 

A 


q,s 


q,s 


q,s 

(x) 

°P, q, n X ' 


C* (x , u>) 
m,n v ' ’ 


dC-jjAa. 

c 

c 0 


D 

D 


m, n 


6 ij 


F 

F 

f 

t 


m, s 

y 

0 

q 

o 

a,q 


aspect ratio 

Fourier coefficient of axial distortion velocity 
Fourier coefficient of circumferential distortion profile 
number of blades 

Fourier coefficient of radial distortion velocity 
blade height 

Fourier coefficient of circumferential distortion velocity 

defined by eq. (51) 

defined by eq. (18) 
lift slope 

rotor chord length 
speed of sound 

compressibility correction factor defined by eq. (32) 
defined by eq. (16b) 

Fourier coefficient of radial distortion profile 
defined by eq. (39) 
viscous stress tensor 

Fourier coefficient of f a 

Fourier coefficient of f® 

dimensionless circumferential distortion profile 
force per unit area acting on cascade 

force on individual blade in cascade 


G 

g 

H P, q, s 


force per unit projected area acting on cascade 
Green’s function 
radial distortion profile 
defined by eq. (43) 
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H. 


( 0 ) 

p, q, s 


*p 

K 


p, q, S 


m, n 


(k) 


l 

M 

M r 

M t 

N 

P D 

P Q 

P(w) 

*P 

P’ 


**P 

PQ 


defined by eq. (42) 

ratio of maximum distortion velocity to mean axial velocity U 

AT* 

axial component of the sound intensity in the p harmonic 
defined by eq, (44) 
wave number, oj/Cq 

V _ q ~2 

k ~ @ K m n 

Fourier coefficient of l 

lift force per unit length of blade 

axial Mach number, U/Cq 

relative Mach number, U p /cq 

rotor tip Mach number, U^/c q 

positive integer, 1,2,. . . 

dipole contribution to sound pressure spectrum 

quadrupole contribution to sound pressure spectrum 

sound pressure spectrum 

J.U 

total acoustic power radiated upstream in the p harmonic 
pressure measure above ambient pressure 

J-U 

p harmonic of acoustic pressure fluctuation 
ambient pressure 


n defined by eq. (16a) 

R mean radius 

S a , S b vertical duct surface area 

S f rotor surface area (relative to fixed coordinator) 

S® rotor surface (relative to rotor-locked coordinates) 

S(Oq) Sears function 

T large time interval (taken equal to infinity at end of analysis) 

tJ. LighthilFs stress tensor based on relative velocity 

t time (associated with observation point) 

U axial flow velocity 
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u * 

u* 


U. 


v i 


1 


w. 


X 

x 

y 

z 

z 


a i 


0 

h 


y 

A 


A 


0 


6 

5.. 

i] 

6(x) 


V 

°p, q, s 

9 


K m, n 


X 


s,q 


relative velocity, - ^tf 2 + U 2 
cascade (rotor) velocity 
distortion velocity 

phase-locked rotor velocity field relative to rotor blades, i = 1,2 
total fluid velocity 
relative velocity, v^ - 6 U U 
p harmonic of acoustic velocity 
phase-locked rotor velocity field 
coordinates fixed to rotor 
coordinates of observation point 
coordinates of source point 
complex variable defined in terms 
complex variable defined in terms 
defined by eqs. (17) 
y 1 - M 2 
■^1 - M 2 

V 1 -"? 

circulation about blade 
stagger angle 
interblade spacing 
defined by eq. (28) 
circumferential distance, BA 
Krone eke r delta 
delta function 

rotor- locked coordinate system (eq. (12)) 
defined in eq. (64) 
work coefficient, -AWg/U^ 
eigenfunctions of wave operator 
defined by eq. (69) 


of X coordinates by eq. (26) 
of rf coordinates by eq. (31) 
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u(r) region of duct exterior to blades 

i>q v(t) at a fixed instant of time 

Pq background density 

p' density fluctuation in acoustic field 

a\ cr" half- width of Gaussian/step profile 

o reduced frequency given by eq. (57) 

M. 

a dimensionless radial hieght of flow distortion 

r time (emission) 

<i> n cross- duct eigenfunctions 
i// defined by eq. (65) 

O shaft rotational frequency, 2nV^/5 

a) frequency 

Subscripts: 

i,j,k integer 1,2,3 

n integer 0, 1,2, . . . 

s integer 0, 1,2, . 

p, q integer 0, ±1, ±2, . . . 

a integer 1,2,. . . 

Superscripts: 

± upstream propagation/downstream propagation 

t relative to coordinates moving with axial velocity U 

C~ ) complex conjugate 
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TABLE I. - DESIGN PARAMETERS AND SUBSONIC OPERATING CONDITIONS 


FOR 51-CENTIMETER (20- IN. ) MODEL FAN (ROTOR 1 - MODEL 1) AT 
5 PERCENT OF SPAN FROM TIP 

[interblade spacing, A, 3. 52 cm (1. 385 in. ); ratio of blade length to interblade 
spacing, b/A, 3.9; number of blades, B, 45; solidity, c/A, 1.3.] 


Corrected fan speed, 
percent of design 

Axial Mach 
number, 

M 

Relative Mach 
number, 

M r 

Work coefficient, 

e 

80 

0.525 

0.973 

0.129 


.463 

.938 

.307 


.382 

.896 

.322 

70 

0.480 

0. 864 

0. 0842 


.40 

.818 

.302 


.31 

.771 

.327 

60 

0.446 

0.759 

0. 0508 


.340 

.696 

.294 


.268 

.661 

.331 


TABLE II. - DESIGN AND OPERATING PARAMETERS FOR FULL-SCALE FAN STAGES AT 5 PERCENT OF SPAN 

FROM TIP (OPERATED WITH DESIGN NOZZLE) 


Fan 

Number of 
blades, 

B 

Blade length, 
b 

Interblade spacing, 
A 

Solidity, 

c/A 

Blade Mach 
number, 

M t 

Axial Mach 
number , 

M 

Corrected fan 
speed, 
percent of 
design 

Work coeffi- 
cient, 

6 

cm 

in. 

cm 

in. 

QF-1 

53 

49 

19.35 

10.6 

4. 15 

1.34 

0.86 

0. 45 

90 

0.46 

QF-3 

53 

49 

19.35 

^ 6 

4.15 

1.34 

.86 

.45 

90 

.375 

A 

40 

49.7 

19.62 

3 9 

5.47 

1.45 

.815 

.40 

80 

.425 

B 

26 

49.7 

19.62 

21.5 

8.45 

1.30 

.814 

.40 

80 

.425 
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TABLE m. - COMPARISON OF THEORY WITH FULL-SCALE FAN 
STAGES (WITH OPERATING CONDITIONS OF TABLE H) 


Fan 

Experimental 

Theory 


Power radiated in fundamental tone, dB 

QF-l 

147 

147 

QF-3 

146 

145 

A 

147 

142 

B 

149 

141 
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(a) Annular duct. 
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(b) Annular duct unrolled - ractanguiar duct. 
Figure 2. - Duct cross section. 



Figure 3. - Configuration of fan in rectangular geometry. 




25, 



Figure 6. - Effect of blade number on power in fundamental tone. Half 
Half-width of Gaussian profile, o’, 0. 1; dimensionless width of step 
profile, a, 0.2; axial Mach number, M = 1/2 M t ; solidity, c/A, 1.3; 
work coefficient, 8, 0.4 


Ha If -width of 
Gaussian 
profile, 
o' 



Figure 7. - Effect of distortion size on tone power level. Dimension- 
less width of step profile, o, 0.2; Mach number, M = 1/2 M t ; sol- 
idity, c/A, 1.3; number of blades, B, 40; work coefficient, 0, 0.4 
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Figure 8. - Front -quarter view of rotor 1 - model 1. 
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Fundamental tone, sound power level in 50-Hz filter, dB re 10” 1:> W 


Corrected fan 
speed, 



Figure 9. - Comparison of radiated sound power with 51-centimeter (20-in. ) mod 
el fan (rotor 1 - model 1). Width of Gaussian profile, 2o', 0.02; ratio of maxi- 
mum distortion velocity to mean axial velocity, I, 0.015. 



Figure 10. - Full-scale fan test installation. 
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